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Abstract 

Symmetry properties of the nonlinear heat conductivity equations of the general form u t — [E(x, u)u x ] -\ 
H{x, u) is studied. The point symmetry analysis of these equations is considered as well as an equivalence 
classification which admits an extension by one dimension of the principal Lie algebra of the equation. 
The invariant solutions of equivalence transformations and classification of the nonlinear heat conductivity 
equations among with additional operators are also given. 



1 Introduction 



o 

Q 
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C3 i Symmetry properties of mathematical models of heat conductivity and diffusion processes [8] are tradi- 
£h ! tionally formulated in terms of nonlinear differential equations which often envisage us with difficulties in 
studying. To solve this problem, symmetry methods play a key role for finding their exact solutions, similar 
(N solutions and invariants [2l [4l [TBI [T9l 120] . 

, In this study, we generalize the study of a class of the nonlinear heat conductivity equations (HCEs) 
which has been recently studied in some special cases [13 ^ 119 ^ 120], We deal with the class of nonlinear heat 
t" — ■ conductivity equations of the general form 

cn ; 

OV HCE: u t = [E(x,u)u x ] x + H{x,u), (1.1) 

o ; 

in which we assumed that E, H are sufficiently smooth functions, E x , E u , H x , H u ^ 0, u is treated as the 
• • dimensionless temperature, t and x are the dimensionless time and space variables and E is the thermal 
. ^ conductivity. 

^ ■ The Lie point symmetry in linear and nonlinear special cases of our problem was investigated. For 
instance the problem for the case in which E just depends to u and E u = corresponding to the linear 
case, the case where the class of nonlinear one-dimensional diffusion equations when E = E(u), H = 0, the 
class of diffusion-reaction equations when E = E(u),H = const, the case which the thermal conductivity 
is a power function of the temperature and additional equivalence transformations, conditional equivalence 
groups and nonclassical symmetries have all investigated and listed in Table 1 of [20]. The point symmetry 
group of nonlinear fin equations of class (jl.lB were considered in a number of papers. Moreover in |13j . 
preliminary group classification of nonlinear fin equation was studied in the general form. 

Class of Eqs. (jl.ip generalizes a great number of the known nonlinear second order equations describing 
various processes in biology, ecology physics and chemistry (see [2] and also [3] with references therein). 
Nonlinear heat equations in one or higher dimensions are also studied in literature by using both symmetry 
as well as other methods [5j [6] (an account of some interesting cases is given by Polyanin p2]). But for 
the first time, we generalize the equations described above to the nonlinear heat conductivity equations in 
the form (|l.lj) to investigate their symmetry properties. In order to determine more symmetry of HCEs, 
after finding the point symmetry group, we use preliminary group classification to find different cases of 
one-dimensional extension of the symmetry algebra. 

The more general class of HCEs is the nonlinear heat conductivity equations of the form 



u t = F(t, x, u, u x )u xx + G(t, x, u, u x ), 



(1.2) 



which admits non-trivial symmetry group. The group classification of (jl.2p is presented in some references 
[UHI?]. However, since the equivalence group of (|1.2j) is essentially wider than those for particular cases, the 
results of [U [12] cannot be directly used for symmetry classification of particular ones. Nevertheless, these 
results are useful for finding additional equivalence transformations in the class of our problem. Therefore in 
contrast to the above works, in the last two sections of this paper, we study group classification of Eq. (jl.ip 
under equivalence transformations in the general case. Furthermore, a number of nonlinear invariant models 
which have nontrivial invariance algebras are obtain. 

From [SI] we know that if the partial differential equation possesses non-trivial symmetry, then it is 
invariant under some finite-dimensional Lie algebra of differential operators which is completely determined 
by its structural constants. In the event that the maximal algebra of invariance is infinite-dimensional, 
then it contains, as a rule, some finite-dimensional Lie algebra. Also, if there are local non-singular changes 
of variables which transform a given differential equation into another, then the finite-dimensional Lie 
algebra of invariance of these equations are isomorphic, and in the group-theoretic analysis of differential 
equations such equations are considered to be equivalent. To realize the group classification, we use of the 
proposed approach consists in the implementation of an algorithm explained and performed in references 
[H El [15| 118] , For this goal, our method is similar to the way of for the nonlinear wave equation 
utt = f(x, u)u xx + g(x, u). 

In the next section, we concern with the problem of finding point symmetry group of Eq. (jl.ip . In 
reminded sections of the paper we find some further symmetry properties of Eq. (jl.ip by use of equivalence 
transformations and an extension by one dimension of the principal Lie algebra of the equation. 

2 Lie point symmetries 

In this section, our study is based on the method of |14| for Lie infinitesimal method, we are concerning 
with group classification of HCEs by the point transformations group. 

An equation of class (jl.ip is a relation among with the variables of 2-jet space J 2 (R 2 ,R) with (local) 
coordinate 

(t,x,u,u t ,u x ,uu,u tx ,u xx ), (2.3) 

where this coordinate involving independent variables t, x and dependent variable u and derivatives of u in 
respect to t and x up to order 2 (each index will indicate the derivation with respect to it, unless we specially 
state otherwise). Let M. be the total space of independent and dependent variables resp. t,x and u. The 
solution space of Eq. (jl.ip . (if it exists) is a subvariety Sa C J 2 (R 2 ,R) of the second order jet bundle of 
2-dimensional sub-manifolds of A4. Point symmetry group on Ai is introduced by transformations in the 
form of 

t = Q(t, x,u), x = H(t, x, u), u = Q(t, x, u), (2.4) 

for arbitrary smooth functions (p,x^- Also assume that the general form of infinitesimal generators is 

Odd 
Y := x, u) — + T(t, x, u) — + cp(t, x, u) — , (2.5) 

when coefficients are arbitrary smooth functions. These infinitesimals signify the Lie algebra C of the point 
symmetry group G of Eq. (jl.ip . The second order prolongation of X [14j [T5] as a vector field on J 2 (R 2 , R) 
is as follows 

yW> : = v + (p td +ipX d +iptt d + ^ d + 8 

OUt OU x OUu OUt x OU xx 

where rf,rf and rf l , rf x , rj xx are arbitrary smooth functions depend to variables t, x, u, ut,u x and (|2.3p resp. 
These coefficients are introduced as following 

7] J = Vj{Q) + iuj, t + TUj )X (2.7) 



where T> is total derivative, J is a multi-index with length 1 < | J\ < 2 of variables t, x and Q = u — ^ut — ru x 
is characteristic of v [14] . According to |14j . u is a point infinitesimal generator of Eq. if and only if 
Y^fEq. (jl.ip ] = when Eq. (jl.ip is hold. By applying on the equation we have the following equation 

t[E x u xx + E u u x u xx + E xx u 2 x + E uu u x + E xx u x + E ux u 2 x + H u u x + H x ] + (p[E u u xx + E uu u 2 x 
+E UX u x + H u ] -(p t + (f x [2E u u x + E x ] + E <p xx = 0, whenever Eq. (TTT]) is satisfied. (2.8) 

In the extended form of the latter equation when we consider ut = [E(x, u)u x ] x + H (x, u), functions £, r 
and tp only depend to t,x,u rather than other variables ut,u x ,utt,ut x and u xx , hence the equation will be 
satisfied if and only if the individual coefficients of the powers of ut,u x and their multiplications vanish. 
This tends to the following over-determined system of determining equations 

E£ x = 0, E£ u = 0, E u t + 2Et u = 0, 

E t uu + E u t u - E uu t = 0, E x t + E u <p - 2Et x + E£ t = 0, 

H x (t + ip) - (p t + E x ip x + E ip xx + H(£ t - <fu) = 0, 

E u £ t + E(if uu - 2t ux ) - 2E U t x + E u (p u + E uu if + 2E UX r = 0, 

(E xx + 2H u )t + E ux ip + T t + 2E U <p x + E x {& - t x ) + E{2 ip ux - r xx ). (2.9) 
The general solution to differential equations (12.9P for £, r and </3 is 

^(t,x,n)=c, T(t,x,n) = 0, y?(t, x, u) = 0, (2-10) 

for arbitrary constant c. 

Theorem 1. A complete set of all infinitesimal generators of the HCE ([QP up to point transformations 
admits the structure of one- dimensional Lie algebra C = (Jj)- 

It is well-known that the existence of a non-fiber-preserving symmetry usually indicates that one can 
significantly simplify the equation by some kind of hodograph-like transformation interchanging the inde- 
pendent and dependent variables. Since the kernel of maximal Lie algebra of the HCE (jl.ip is C = (m), 
thus we can not use this advantage for simplifying HCEs. 

Furthermore, according to the statements of page 209 of |14j . we conclude that 

Corollary 2. A system of HCEs of class can not be reduced into an inhomogeneous form of a linear 

system. 



3 Equivalence transformations 

In this section, we follow the method of Ovsiannikov [15j for partial differential equations. His approach is 
based on the concept of an equivalence group, which is a Lie transformation group acting in the extended 
space of independent variables, functions and their derivatives, and preserving the class of partial differential 
equations under study. It is possible to modify Lie's algorithm in order to make it applicable for the 
computation of this group [H [15l [18] . Next we construct the optimal system of subgroups of the equivalence 
group. 

An equivalence transformation is a non-degenerate change of the variables t,x,u taking any equation of 
the form (jl.ip into an equation of the same form, generally speaking, with different E(x, u) and H (x, u). The 
set of all equivalence transformations forms an equivalence group G. We shall find a continuous subgroup 
Gc of it making use of the infinitesimal method. 

We investigate for an operator of the group Gc in the general form 

Odd d d 

Y := x, u) — + r(t, x,u) — + <p(t, x,u) — + x(t, x, u, E > h )-g^ + v(t, x, u, E,h)—. (3.11) 



from the invariance conditions of Eq. (jl.ip written as the system 

u t = [E(x,u)u x ] x + H(x,u), E t = H t = 0, (3.12) 

where we assumed that u, E, H are differential variables: u on the base space (t, x) and E, H on the total 
space (t,x,u). Also in Eq. (|3.1ip the coefficients are dependent to t,x,u and the two last ones, in addition, 
depend to E, h. The invariance conditions of the system (|3. 12[) are 

Y [u t - [E(x, u)u x ] x - H(x, u)) =0, Y [E t ] = Y [H t ] = 0, (3.13) 

where 

t 9 _ d TT d f d T d .. d t d . „ jS 

Y:=Y + ^^T t + ^ 8^ + V 0^ x +X 0E t +X dE x + * 8E u +ri 8H t - (3 ' 14) 

is the prolongation of the operator (13. lip . Coefficients rj J for multi-index J (with length 1 < | J\ < 2) have 
given in section 2 and by applying the prolongation procedure to differential variables E, H with independent 
variables (t, x, u) we have 



(3.15) 



X 1 = Pj(x) - E t Vj(0 - E x Vi(t) - E u VM) = Vj( x ) - E x V T (r) - E u Vj(^), 
rf = V t ( V ) - H t V t (0 - H x V t (ri) - H u V t (4>) = V t (rj) - H x V t (r) - H u V t (<j>), 
where I varies on variables t,x,u and 

~ d d d 

Vl --=dI +El dE +Hl dH- 
Substituting (|3.14p in (|3.13|) we tend to the following system 

X u xx +r ] -i P t + [2E U u x + E x ] ip x + E ip xx + x x u x + X u u 2 x = 0, (3.16) 

X I = 0, rf = 0. (3.17) 

Replacing relations r] J (for multi-index J with length 1 < |J| < 2) and X^iX i n Eqs. (13.16p - (|3.17p and 
then introducing the relation ut = [E(x,u)u x ] x + H(x,u) to eliminate ut, we have three relations which are 
called determining equations. The two last ones are the determining equations associated with Eqs. (|3.17p . 
i.e., 

Xt - E x T t - E u (pt = 0, rj t - H x r t - H u ip t = 0. (3.18) 

But these relations must hold for arbitrary variables E x , E u , H x , H u of the jet space and this fact results in 
the following conditions 

T t = 0, <pt = 0, Xt = 0, r]t = 0, (3.19) 

so, we find that 

£ = £,{t,x,u), r] = rj(x,u), tp = ip(x,u), x = x(x,u,E,H), rj = rj(x, u, E, H). (3.20) 

But adding these conditions to the first determining equation, knowing that ut,u x ,uu,ut x ,u xx are considered 
to be independent variables, we lead to the following system of equations 

X + E[£ t -2T X ] =0, Eif uu + X u = 0, £t-2T x + X E = 0, 

E [2 ip ux - t xx ] + Xx = 0, ip x = 0, xh = 0, (3.21) 

ri + H[tt-<Pu] + E<p xx = 0, E£ x = 0, E£ u = 0, Et u = 0. 

This system follows 

^ = (2C1 + C 2 ) t + C 3 , T = ClX + C4, ip = c 5 u + c 6 , 

(3.22) 

X = c 2 E, r] = (c 5 - 2 ci - c 2 ) H, 



A 



for arbitrary constants c\,---,cq. Therefore the class of Eqs. (jl.ip has an infinite continuous group of 
equivalence transformations generated by infinitesimal operators 

Y ® Y ® Y ® Y 2t ^ ^ 2H ® 

1 dt : 2 dx' 3 du 4 dt dx OH' (qoi\ 

8^8,8 „ 8,8 [6 '"> 

Y * = - t M + E aE + H W Y '- , 'Tu + H W 

Moreover, in the group of equivalence transformations are included also discrete transformations, i.e., re- 
flections 

1 1 — > -t, u i — y-u, E\ — ► —E, H \ — ► -H. (3.24) 

The communication relations between these vector fields is given in Table 1. The Lie algebra C := (Yj : 
i = 1, ■ ■ ■ , 6 ) is solvable since the descending sequence of derived subalgebras of C: C D C^ 1 ' = { Yj_, Y<i, Y3 ) D 
£( 2 ) = {0}, terminates with a null ideal. But for each v = Y^i v i¥i an d w = Ylj w jYj m A its Killing form: 

K(v, w) = tr(ad(u) o ad(w)) = 5 04 64 — 2 (04 65 + 05 64) + 05 65 + a% (3.25) 

is degenerate. Hence C is neither semisimple nor simple. 



Theorem 3. Let Gi be the one-parameter group (flow) generated by Yi, then we have 

G\ : (t, x, u, E, h) 1 — > (t + s, x, u, E, H), G2 ■ (t, x, u, E, h) 1 — > (t, x + s,u, E, H), 

G 3 : (t, x, u, E, h) 1 — ► (t, x,u + s, E, H), G 4 : (t, x, u, E, h) 1 — > (t e 2s ,x e s ,u, E, H e~ 2s ), (3.26) 

G5 : (t,x,u,E,h) 1 — > (te~ s ,x,u,Ee s ,H e s ), Gq : (t,x,u, E,h) 1 — ► (t,x,ue s , E, H e s ), 

when s is an arbitrary parameter. Moreover, if u = f(t,x) for functions E and H be a solution of the HCE 
so are 

ui = f(t + s,x), u 2 = f(t,x + s), u 3 = f(t,x,u) - s (3.27) 

for the same functions E and H, U4 = f(te 2s ,xe s ) for the same functions E and H = H e~ 2s , U5 = 
fit e~ s , x) for E = E e s and H = H e s and also uq = e~ s fit, x) for the same E and H = H e s . 

4 Preliminary group classification 

In many applications of group analysis, most of extensions of the principal Lie algebra admitted by the 
equation under consideration are taken from the equivalence algebra These extensions are called £— 
extensions of the principal Lie algebra. The classification of all nonequivalent equations (with respect to 
a given equivalence group Gg) admitting ^-extensions of the principal Lie algebra is called a preliminary 
group classification [9] . We consider the algebra £ spanned on operators (|3,23p and use it for a preliminary 
group classification. 

It is well-known that the problem of classifying invariant solutions is equivalent to the problem of clas- 
sifying subgroups of the full symmetry group under conjugation in which itself is equivalent to determining 
all conjugate subalgebras [HI [15]. The latter problem, tends to determine a list (that is called an optimal 
system) of conjugacy inequivalent subalgebras with the property that any other subalgebra is equivalent to 
a unique member of the list under some element of the adjoint representation i.e. Ch Ad(<?) Ch for some g 
of a considered Lie group. Thus we will deal with the construction of the optimal system of subalgebras of 
C. 

The adjoint action is given by the Lie series 

Ad(exp(s Yi))Yj =Y j -s [Y^Yj] + ^ [Y h [5* Y-]] - • • • , (4.28) 

where s is a parameter and i,j = ,6. The adjoint representations of £ is listed in Tables 2; it consists 
the separate adjoint actions of each element of C on all other elements. 



Table 1: Commutator table 
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Table 2: Adjoint table 
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Theorem 4. An optimal system of one- dimensional Lie subalgebras of the nonlinear HCE M-l\) is provided 
by those generated by 

1) A 1 =Y 1 = d t , 

2) A 2 = Y 2 = d x , 



u • 



3) A 6 = Y 3 = 3, 

4) A 4 = Y A = 2td t + d x -2Hd H , 

5) a 5 = y 5 = -td t + Ed E + Hd H , 

6) a 6 = y 6 = ud u + Hd H , 

7) A 7 = ±Y X + Y 2 = ±d t + d x , 



A 8 = ±Y X + Y 3 = ±dt + a, 



•0.1 



9) A 9 = ±Yi + Y 6 = ±dt + ud u + Hd H , 



10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 



A 10 =Y 2 + Y 3 = d x + d u , 

A 11 = ±Y 2 + Y A = 2td t ± X + d x - 2Hd H , 

A 12 = ±Y 2 +Y 5 = -td t ±d x + Ed E + Hd H , 

A 13 =Y 2 + Y 6 = d x +ud u + Hd H , 

A M = ±Y 3 + Yi = 2td t + d x ±d u - 2Hd H , 

A 15 = ±Y 3 + Y 5 = -td t ±d u + Ed E + Hd H , (4.29) 

A 16 = ai Y 4 + Y 5 = (2ai - l)t d t + a 1 d x + Ed E - (2a x - l)H 8 H , 

A 17 = a 2 Yi + y 6 = 2 a 2 t d t + a 2 d x + ud u - (2a 2 - l)H %, 

A 1S = f3 1 Y 5 + Y 6 = -^t 8 t + ud u + f3 1 Ed E + {fix + 1)H d H , 

A 19 = ±Yx +Y 2 + Y 3 = ±dt + d x + d u , 

A 20 = ±Y X +Y 2 + Y 6 = ±dt + d x + ud u + H8 H , 

A 21 = ±Y 2 ± Y 3 + Yi = 2t 8 t ± d x + d x ± 8 U - 2H8 H , 

A 22 = ±Y 2 ±Y 3 + Y 5 = -tO t ± 8 X ±d u + Ed E + Hd H , 



±Y 2 + a 3 Y 4 + Y 5 = (2a 3 - l)i 8 t + (a 3 ± + E d E — (2a 3 - 1)# 
F 2 + a 4 r 4 + ie = 2a 4 * ft + (a* + l)9a, + u d u - (2a 4 - 1)# d H , 
±Y 2 + (3 2 Y 5 + Y 6 = -p 2 td t ±d x + ud u + f3 2 E8 E + (/3 2 + l)Hd H , 
±Y 3 + a 5 Y A + Y 5 = (2a 5 - l)t d t + a 5 d x ± d u + E d E - (2a 5 -l)Hd H , 
a 6 n + h Y b + Y & = (2a 6 - fa)t 8 t + a & d x + ud u + fi 3 E 8 E - (2a 6 -fa-l)Hd H , 
±Y 2 ±Y 3 + a 7 Y 4 + Y 5 = (2a 7 - l)td t + (a 7 ± 1) ± <9 U + £ - (2a 7 - 
±y 2 + a 8 n + /3 4 F 5 + Y 6 = (2a 8 - /3 4 )t 8 t + (a 8 ± 1) d x + ud u + p 4 E3 E - (2a 8 - /3 4 - H , 

/or nonzero constants cti,(3j (1 < i < 8, 1 < j < 4J. 

Proof. Let /I is the symmetry algebra of Eq. (|l.ip with adjoint representation determined in Table 2 and 

Y = a x Yx + a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + a 5 Y 5 + a 6 Y 6 (4.30) 

is a nonzero vector field of C. We will simplify as many of the coefficients ctj as possible through proper 
adjoint applications on Y. We follow our aim in the below easy cases. 

Case 1 At first, assume that 7^ 0. Scaling Y if necessary, we can consider to be 1 and so follow the 
problem with 

Y = ax Yx + a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + a 5 Y 5 + Y 6 . (4.31) 

Case la According to Table 2 in the case which a 4 7^ 0, if we act on Y by Ad(exp(^- Y 2 )), the coefficient 
of Yx can be vanished: 

Y' = a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + a 5 Y 5 + Y 6 . (4.32) 

Then for 03 ^ we apply Ad(exp(a 3 Y 3 )) on Y' to cancel the coefficient of Y 3 (it is automatically hold 
for a 3 = 0): 

Y" = a 2 Y 2 + a 4 Y 4 + a 5 Y 5 + Y 6 . (4.33) 
Case la-1 If in addition a 2 7^ 0, we can act Ad(exp(±ln(— ) I4)) on y" to change a 2 to ±1: 

F w = ±y 2 + a4i4 + a5^5 + >6- (4.34) 

Then for the case which 05 7^ 0, we can not simplify 1"'" any more. This introduce part 29 of the 
theorem. 

Also when 05 = we tend to part 24 of the theorem. 

Case la-2 In Case la, let we consider a 4 of Y" is equal to zero. Then for 

Y "" = a 4 Y 4 + a 5 Y 5 + Y 6 , (4.35) 

where 05 7^ more simplification is impossible and this results in part 27, and where 05 = it results 
in part 17 of the theorem. 

Case lb Let the coefficient a 4 in Y is zero. So we have the following new form of Y 

Y = axYx + a 2 Y 2 + a 3 Y 3 + a 5 Y 5 + Y 6 . (4.36) 

Case lb-1 In the case which 05 7^ 0, by applying Ad(exp(— ^7) Yx)) on Y, the coefficient a 4 will be vanished: 

Y = a 2 Y 2 + a 3 Y 3 + a 5 Y 5 + Y 6 . (4.37) 

Furthermore, if we act Ad(exp(^|) Y 3 )) on the latter form, so the coefficient a 3 will be zero. When 
a 2 7^ 0, applying Ad(exp(± In — ) Y 4 )) any one-dimensional subalgebra generated by Y is equivalent 
to one generated by ±Y 2 + 05 Y5 + Yq which introduce part 25 of the theorem for constant 7^ 0. For 
the other case which a 2 = 0, part 18 are given. 



23) A 2i 

24) A 24 

25) A 25 

26) A 26 

27) A 27 

28) A 28 

29) A 29 



Case lb-2 In Y, let as = 0. In this case, by applying Ad(exp(2|) Yq)) we can make the coefficient of Y 3 
equal to zero: 

Y = a 1 Y 1 + a 2 Y 2 + Y 6 . (4.38) 

Case lb- 2-1 For a 2 + the action of Ad(exp(± In ^) Y 4 )) on F results in the form ±^Y 1 + Y 2 + Y 6 . 
Then if a\ 7^ 0, we can apply Ad(exp(± In ^2-) Y 5 )) on it to have the form ±Y]_ + Y 2 + Yq that is part 
20 of the theorem. If we consider the case in which a\ = 0, then we tend to part 13. 

Case lb- 2-2 In Y assume that a 2 = 0. In this case the condition a-y 7^ and the action of Ad(exp(± In ^-) Y" 5 )) 
show that the simplest form of Y is similar to part 9, while the condition a\ = leads to part 6 of the 
theorem. 

Case 2 The remaining one-dimensional subalgebras are spanned by vector fields of the form Y with a§ = 0. 
Case 2a If as 7^ then by scaling Y, we can assume that 05 = 1: 

Y = a 1 Y 1 + a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + Y 5 . (4.39) 

Case 2a-l Now by the action of Ad(exp(^- Y 2 )) on Y where 04 / 0, we can cancel the coefficient of Y\\ 

Y' = a 2 Y 2 + a 3 Y 3 + a 4 Y 4 + Y 5 . (4.40) 

Case 2a-l-l Then for 03 / by applying Ad(exp(± In ^ Yq)) on Y' the coefficient of Y 3 can be ±1. Thus 
if a 2 7^ the action of Ad(exp(± In ^ Y4)) recommends part 28 of the theorem. Moreover, where 
a 2 = we tend to part 26. 

Case 2a-l-2 In case 2a-l, let a 3 = 0. So for a 2 / applying Ad(exp(± In ^ Y 4 )) on Y' one can lead to 
part 23 of the theorem, whereas for a 2 = 0, part 16 are achieved. 

Case 2a-2 Suppose that a 4 in Y is equal to zero, then for a\ 7^ by acting Ad(exp(— 21 Y 5 )) and also for 
di = we lead to the following form 

Y" = a 2 Y 2 + a 3 Y 3 + Y 5 . (4.41) 

Case 2a-2-l For 03 / we can act Ad(exp(± In Yq)) on Y" to change the coefficient of Y3 to be equal 
to either +1 or —1. Hence for a 2 7^ by applying Ad(exp(±ln Y4)) and for a 2 = we find part 22 
and part 15 resp. 

Case 2a-2-2 Assume that in Y, a 3 is zero, then for a 2 7^ by applying Ad(exp(±ln ^ Y 4 )) and for a 2 = 
we tend to part 12 and part 5 of the theorem resp. 

Case 2b In the situation of Case 2, suppose that as = 0. 

Case 2b- 1 If in addition we assume that 04 7^ 0, then if necessary we can let it equal to 1: 

Y = a 1 Y l + a 2 Y 2 + a 3 Y 3 + Y 4 . (4.42) 
So for a\ 7^ by applying Ad(exp(|^ Y 2 )) and also for a\ = 0, we have the following form for Y 

Y' = a 2 Y 2 + a 3 Y 3 + Y 4 . (4.43) 



O 



Case 2b-l-l If 03 7^ 0, the action of Ad(± lnexp(^- Y$)) shows that we can make the coefficient of I3 equal 
to ±1. Then for 02 7^ by applying Ad(exp(±ln — Y4)) and for a 2 = resp. part 21 and part 14 of 
the theorem are given. 



Case 2b-l-2 In Y' , suppose that a 3 = 0. The simplest possible form of Y is equal to part 11 after taking 
a 2 7^ and after acting Ad(±lnexp(— Y4)). Moreover, when 02 = the simplest case is equal to part 
4. 

Case 2b-2 With conditions of Case 2b, in addition let 0,4 = 0: 

Y = a 1 Y 1 + a 2 Y 2 + a 3 Y 3 . (4.44) 

Case 2b- 2-1 Consider a 3 7^ 0, then by scaling we can make the coefficient of Y equal to 1. By assuming 
a 2 7^ and applying Ad(exp(± In — Y4)) we find the following form 

Y' = ±—Y 1 + Y 2 + Y 3 . (4.45) 

When a\ 7^ we can reduce it to part 19 after acting Ad(— exp(^- Y5)), whereas for a\ = we 
find part 10. If we change our assumption on a 2 and consider it equal to 0, for a\ 7^ by applying 
Ad(± lnexp(^j- Ye,)) and for a\ = resp. we lead to part 8 and part 3. 

Case 2b- 2-2 Finally if in the same conditions with Case 2b-2, we assume that 03 be zero, by scaling we 
can make the coefficient of Y 2 equal to 1 if a 2 7^ 0: 

Y" = a 1 Y 1 +Y 2 . (4.46) 

If a\ 7^ we can act Ad(±lnexp(— In Y§)) on Y". Then no further simplification is possible and 
then Y is reduced to part 7, while a\ = suggests part 2. The last case occurs when we change the 
condition on a\ to be equal to zero, which recommends part 1 of the theorem. 

There is not any more possible case for studying and the proof is complete (Note that in the group of 
equivalence transformations there are included also discrete transformations, Eq. (|3.24p ). A 

The coefficients E, H of Eq. (jl.ip depend on the variables x, u. Therefore, we take their optimal system's 
projections on the space (x,u,E,h). The nonzero in (x,u)— axis projections of ()4.29p are 

1) Z 1 = A 2 = A 7 = d x , 

2) Z 2 = A 3 = A 8 = d u , 

3) Z 3 = A 4 = d x -2Hd H , 

4) Z* = A 5 = Ed E + Hd H , 

5) Z 5 = A 6 = A 9 = ud u + Hd H , 

6) Z 6 = A 10 = A 19 = d x + d u , 

7) Z 7 = A 11 =±d x + d x -2Hd H , 

8) Z 8 = A 12 = ±d x + Ed E + Hd H , 

9) Z 9 = A 13 = A 20 = d x + ud u + Hd H , 

10) Z 1 

11) Z 1 

12) Z lz = A w = aid x + Ed E - (2ai - l)H d H , (4.47) 

13) Z 1 

14) Z 1 

15) Z 1 

16) Z 1 



10 


= A 14 


= d x ±d u -2Hd H , 


11 


= A 15 


= ±d u + Ed E + Hd H , 


12 


= A 16 


= a 1 d x + Ed E -(2a l -l)Hd H , 


13 


= A 17 


= a 2 d x + ud u - (2a 2 - 1)H d H , 


14 


= A 18 


= ud u + p 1 Ed E + ((3 1 + l)Hd H , 


15 


= A 21 


= ±d x + d x ±d u -2Hd H , 


16 


= A 22 


= ±d x ±d u + Ed E + Hd H , 



17) Z 17 = A 6 = (a 3 ± 1)3, + £<9 £ - (2a 3 - l)i2"0jj, 

18) Z 18 = A 24 = (04 + l)a K +«5 u - (2a 4 - 1)£T Sh, 

19) Z 19 = A 25 = ±d x + ud u + h Ed E + (/3 2 + 1)H d H , 

20) Z 20 = A 26 = a 5 8 x ±d u + Ed E - (2a 5 - 1)# d H , 

21) Z 21 = A 27 = a 6 5 :c + na u + /33 J E3 i; -(2a 6 -/33-l) J ffa // , 

22) Z 22 = A 2S = (a 7 ±l)d x ±d u + Ed E - (2a 7 - 1)H d H , 

23) Z 23 = A 29 = (a 8 ± 1) d x + « d u + /3 4 E - (2a 8 - /3 4 - 1)# d H , 

From paper 7 of [S] we conclude that 

Proposition 5. Let C m := (Yi : i = l,---,m) be an m~dimensional algebra. Denote by A 1 (i = 
1, < s < m, s € N) an optimal system of one-dimensional subalgebras of C m and by Z % (i = 

1, ■ • • ,t, < t < s, t 6 N) the projections of A 1 , i.e., Z % = pr(^4 l ). // equations 

F = F(x, u), G = G(x, u), (4.48) 
are invariant with respect to the optimal system Z l then the equation 

u t = {F(x,u)u x ) x + G(x,u)u, (4.49) 
admits the operators X' 1 = projection of A 1 on (t,x,u). 



Proposition 6. Let Eq. an d the equation 

u t = (F(x,u)u x ) x + G(x,u), (4.50) 

be constructed according to Proposition [7] via optimal systems Z % and Z % resp. If the subalgebras spanned 
on the optimal systems Z % and Z l resp. are similar in C m , then Eqs. \4-4^ ana " \4-5ty are equivalent with 
respect to the equivalence group G m generated by C m . 

Now by applying Propositions U] and U] for the optimal system (|4.47p . we want to find all nonequivalent 
equations in the form of Eq. (jl.ip admitting f-extensions of the principal Lie algebra Cg, by one dimension, 
i.e, equations of the form (jl.ip such that they admit, together with the one basic operator of C\, also a 
second operator X^ 2 \ In each case which this extension occurs, we indicate the corresponding coefficients 
E, H and the additional operator X^ 2 \ 

We perform the algorithm passing from operators Z l (i = 1, • • • , 23) to E, H and X^ via the following 
examples. 

Let consider the vector field 

Z 4 = Ed E + H3 H , (4.51) 
then the characteristic equation corresponding to Z 4 is 



dE _ dH 
~E~ ~ ~H~' 

which determines invariants. Invariants can be taken in the following form 



(4.52) 



Ii = x, I 2 = u, h = —. (4.53) 

hi 

In this case there are no invariant equations because the necessary condition for existence of invariant 
solutions (see [15], Section 19.3) is not satisfied, i.e., invariants (|4. 53|) cannot be solved with respect to E 
and H since each two of them can not be an invariant function with respect to the third one. 



Table 3: The result of the classification 



N 


Z 








Invariant A 


Equation 


Additional operator 


1 


Z 1 










u t = 


[$ u x ] x + * 


A. ±JL + Ji 


2 


z 2 








a; 


u t = 




a 4. 9 , a 
au' at " r a« 


3 


z 3 








M 


u t = 


r/r\ 1 1 —2(x-\-^) 
[$u x \ x + e 


at ^ ax 


4 


z 5 








X 


u t = 


r t i In ni _L \T/ 

[<$>u x ] x + e lnu+ * 


L at 


5 


z 6 








u — X 


Ut = 


[<su x \ x + w 


a 1 a 1 1 c? 1 
dx "f" a«' ^ at "T" ax i" a« 


6 


z 7 








u 


111 — 


\y? u x \ x ^ c 




7 


z 8 








u 


11+ = 


fe* ±x u T l T + e* ±:E 


11=1 at T ax 


8 


z 9 








lnw — x 


u t = 




ax T " a« ' ^ at ^ ax ^ u a« 


9 


z u 








X ±u 


7/+ = 


\$uJr + e - 2 ( a '±«) 


ST — A- + 

zt at ^ aa: ^ 


to 


z 11 








X 


= 




1 at ^ a« 


11 


z 12 


>1 + 1/2) 






u 


= 


[p^+si 21 ,, 1 +P ^+( 2 -^ x ) 


(2«i-l)*^+«i^ 


12 


z 12 


>1 = 1/2) 






u 


u t = 


[e °i u x \ x + * 


_a_ 

ax 


13 


z 13 


> 2 ^ 1/2) 






In it — — 


u t = 


r t "1 , ( 1 Win 7/4-\Ef1 




14 


Z 13 


>2 = 1/2) 






lnu — 2x 


u t = 


\ sf\ ~ . 1 1 VT* 

"xjx + * 


2*I±al± 2 "al 


15 


z u 


;a ^ -i) 






X 


It J. — 


U3i(*(u)-Hnu) 1 , Jl+/3i)(*Hn«) 


^7|-"al 


16 


Z 14 


;a = -1) 






X 


Ut = 




t — — 11 — 

L at "a M 


17 


z 


;29 x ± d u - 2Hd H ) 




1 


u t = 


= [*(«) u x ] x + e- 2 (*±") 


C) 1 Q 1 r) a _|_ 9 

zt at +z el ± a^ 




715 

Zj 


:±o u - 2Hd H ) 






X 


u t = 


[$« x ] x + e- 2 (*±") 


9 + a _i_ a 
z 1 at 31 an 


1 Q 


7I6 
Zj 








ix i 7/ 


u t = 


[e* ± "« x ] x + e* ± » 


— f 9 ± d ± ^ 

at 3x du 


20 


Z 17 


> 3 + ±1,1/2) 








Lit 


r *+(a 3 ±i)x 1 , * ^tt x 




21 


Z 17 


>3 = 1/2) 








Ut = 


r *+(i/2±i)x 1 , ^ 


Pi 

dx 


22 


Z 18 


;a4^ -1,1/2) 






lnU a 4 +l 


= 


, T , 2a 4 — 1 
- T(b 7/ 1 J_ Q4 + I 33 


2a4t-§- t + (a 4 + l)£+u£ 


23 


Z 1B 


> 4 = -1) 






X 


11+ = 


[*uJ T + e 3 (* +lnu ) 




24 




> 4 = 1/2) 






In u — ji 


Wt = 




1 9t + 2 9i + " 9« 


25 


Z ' 


^ -1) 






1 1 

mu ± x 


U f = 


[ e ft(*±x) Ux ] x + e (ft+l)(*±x) 


-fa St ± Tfx-+ u £ 


ZD 


Z 


= -1) 






In 7i ± x 


Ut = 


[e* ±x u x ] x + * 


4. d \ d 1 9 
6 at ^ Oa: U 


27 


Z 20 


>5 + 1/2) 






-2- ±u 


111 = 


\e* ±u u x l x + e*- (2 -* )x 


(2«5-l)i|+a 5 ^±| I 


28 


z 20 


>5 = 1/2) 






f ±« 


= 


[e* ±u u x ] x + * 


— =b 2 — 

5a; du 


29 


Z 21 


> 6 -ft^l) 






In u — — 




r eJ 8s(*+a 6 x) u i . e * 6 4 3 x 


(2a 6 -(3 3 )t-§- t +a 6 JL+u£; 


30 


Z 2l 


> 6 - & = 1) 






lntt — — 


7/ + 


r e /3 3 (*+a 8 !E) w 1 4. v]/ 


fi fi Pi 


31 


Z 22 


>7 + ±1,1/2) 










r *±« 1 , _(l-2a 7 )(*±u) 
L 1 -' u 'X\x 1 1 


(2a r -l)t^ + (a 7 ±l)^±|j 


32 


Z 22 


>7 = ±1) 






X 


■ = 


r^*±u.. 1 1 „(i±2)(*±w) 


(±2-l)^±|l 


33 


Z 22 


> 7 = 1/2) 






X + 7/ 
1/2±1 ^ " 


u t = 


[e* ±M U X ] X + * 


\2 ^ L ) dx du 


34 


Z 23 


^as 7^ ±1, 2a 8 - 


04 


^1) 


lnU a 8 ±l 


u t = 


,t, 2ao-/9 4 -l 




35 


Z 23 


> 8 7^ ±l,2a 8 - 


04 


= 1) 


lnU a 8 ±l 


u t = 


[ e 4 ($+ln«) + ^ 


*& + («8±l)&+«& 


36 


Z 23 


[as = ±1, 2a 8 - 


04 


^1) 


X 


u t = 


rg8 4 (*-Hn«) Ux j x _|_ e (2Q8-/34-l)(*-Hn«) 




37 


Z 23 


'a 8 = ±1, 2a 8 - 


04 


= 1) 


X 


u t = 


[ e (±2-l)(*+lnu) u j a$ 





1 1 



Considering Z 23 (for as ± 1 ^ 0, 2ag — /3 4 — 1 ^ 0) we have the below characteristic equation 

dx du dE dH 
a 8 ±l ~~u ~ TkE ~ (2a 8 -/? 4 - 

This equation suggest the following invariants 



(4.54) 



h = In u - (as ± J 2 = In - In «, I 3 = In H + — — ^ x. (4.55) 

"8 ± 1 

From the invariance equations we can write 

h = $(/!), h = (4.56) 

They result in the forms 

£ = exp(/3 4 (lnu + $(A))), g = exp ( tt(A) - 2 " 8 - ^ - 1 A (4.57) 

\ «8 ± 1 / 

for invariant A = \nu — (as ± l)x. 

From Proposition [5] applied to the operator Z 23 (for as ± 1 ^ 0, 2ag — 04 — 1 ^ 0) we obtain the additional 
operator 

X® = (2a 8 - /3 4 )t d t + (a 8 ± 1) d x + u d u . (4.58) 

One can perform the algorithm for other Z l s of (I4.47D similarly. The preliminary group classification of 
nonlinear fin equation (jl.ip admitting an extension £2 of the principal Lie algebra L\ is listed in Table 3. 



5 Conclusion 

Symmetry analysis for equations ut = [E(x, u)u x ] x + H(x, u) rather than previous results on special cases 
of this equation equation [13 [20], is carried out exhaustively. Also, equivalence classification is given of the 
equation admitting an extension by one of the principal Lie algebra of the equation. The paper is one of few 
applications of a new algebraic approach to the problem of group classification: the method of preliminary 
group classification. Derived results are summarized in Table 3. 
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